Adiabatic nonlinear waves with trapped particles: III. Wave dynamics 
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The evolution of adiabatic waves with autoresonant trapped particles is described within the 
Lagrangian model developed in Paper I, under the assumption that the action distribution of these 
particles is conserved, and, in particular, that their number within each wavelength is a fixed 
independent parameter of the problem. One-dimensional nonlinear Langmuir waves with deeply 
trapped electrons are addressed as a paradigmatic example. For a stationary wave, tunneling into 
overcritical plasma is explained from the standpoint of the action conservation theorem. For a 
nonstationary wave, qualitatively different regimes are realized depending on the initial parameter 
S, which is the ratio of the energy flux carried by trapped particles to that carried by passing 
particles. At S < 1/2, a wave is stable and exhibits group velocity splitting. At S > 1/2, the 
trapped-particle modulational instability (TPMI) develops, in contrast with the existing theories 
of the TPMI yet in agreement with the general sideband instability theory. Remarkably, these 
effects are not captured by the nonlinear Schrodinger equation, which is traditionally considered as 
a universal model of wave self-action but misses the trapped-particle oscillation-center inertia. 

PACS numbers: 52.35.-g, 52.35.Mw, 52.25.-b, 45.20.Jj 
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I. INTRODUCTION 

Within the geometrical-optics (GO) approximation, 
adiabatic nonlinear waves in collisionless plasma are de- 
scribed conveniently within the average-Lagrangian for- 
malism originally proposed by Whitham [B H]- In our 
Refs. [2113], further called Papers I and II, this formalism 
was restated to also accommodate effects of autoreso- 
nant particles trapped in wave troughs. Specifically, the 
corresponding Lagrangian density £ and the nonlinear 
dispersion relation (NDR) were derived in Papers I and 
II (see also Ref. [5] and references therein), under the 
assumption that the action distribution of these parti- 
cles is conserved, and, in particular, that their number 
within each wavelength is a fixed independent parameter 
of the problem. Here, the evolution of such waves will be 
studied. 

It is commonly believed that, within the GO approxi- 
mation, the wave evolution can be inferred from the NDR 
alone, in generalization of the well-known linear solution 
[5]. This implies, in particular, that £ is deducible from 
the NDR. Since the latter is equivalent (Paper I) to 



£a=0. 



(1) 



where a is the amplitude [7], the assumption hereby is 
that the integration constant associated with Eq. (fTl) is 
insignificant, so one can take 2, — L Ha da. However, 
this does not apply to waves with autoresonant trapped 
particles. First of all, at small enough a the adiabatic 
Lagrangian ceases to exist, because particles start to es- 
cape the wave potential. (Only (5-shaped distributions of 
trapped particles are allowed at a — >■ 0.) Second, as flows 
from Paper I, £ can contain terms which are indepen- 
dent of a [so they cannot be inferred from Eq. (IT])] and 
yet depend on uj and k, thus affecting the wave dynam- 
ics. This is also understood from the fact that trapped 
particles carry fractions of the wave momentum density 



p — k£,i^ and the energy flux density H = —uj£k ;2 : a 
part of which is determined by the phase velocity rather 
than a. Therefore, examining just the NDR is insuffi- 
cient to predict the evolution of these waves, contrary to 
Refs. [51-fT5]. Instead, the complete Lagrangian must be 
used [T5]. 

Here, we show how effects captured by £ but not by 
the NDR render the wave self-action due to autoresonant 
trapped particles unique among other self-action mech- 
anisms. Specifically, we analyze the evolution of one- 
dimensional (ID) nonlinear Langmuir waves with deeply 
trapped autoresonant electrons as a paradigmatic ex- 
ample, illustrating the qualitative physics that is also 
expected for other distributions (except when dynamic 
trapping and detrapping become essential). For a sta- 
tionary wave, tunneling into overcritical plasma |17H19) is 
explained from the standpoint of the action conservation 
theorem (ACT). For a nonstationary wave, qualitatively 
different regimes are realized depending on the initial pa- 
rameter S, which is the ratio of the energy flux carried 
by trapped particles to that carried by passing particles. 
At iS* < 1/2, a wave is stable and exhibits group velocity 
splitting. At 5' > 1/2, the trapped-particle modulational 
instability (TPMI) develops, in contrast with the exist- 
ing theories of the TPMI f9'-fTl] yet in agreement with 
the general sideband instability (SI) theory [50]. Re- 
markably, these effects are not captured by the nonlinear 
Schrodinger equation (NLSE), which is traditionally con- 
sidered as a universal model of wave self-action [6] but 
misses the inertial of the trapped-particle oscillation cen- 
ters (OC). 

The work is organized as follows. In Sec. |ll) we for- 
mulate our analytical model in general. In Sec. |III[ we 
present the wave Lagrangian in a simple form and de- 
rive the GO equations flowing from it. In Sec. IV we 



study waves in plasmas with parameters varying in space 
and time; in particular, wave tunneling into overcritical 



plasma is discussed. In Sec. [Vj we consider pulse prop- 
agation in homogeneous stationary plasma; specifically, 
the group velocity splitting is derived for S < 1/2, and 
the TPMI rate is calculated for S > 1/2 



VI 



we 



In Sec 

explain our results in the context of a more general stabil- 
ity criterion and compare them with other existing the- 
ories. In Sec. |VII[ we summarize our main conclusions 
Also, supplementary material is given in appendixes. 



II. BASIC MODEL 



electrons have (Paper I) 

Hi*) = £■(*) - meuV2, 



(5) 



At) 



where £e is the particle energy in the frame where the 
wave is stationary, u = ui/k is the phase velocity, and nie 
is the electron mass. Here, we will assume, for simplicity, 
that these particles initially reside at (and, due autores- 
onance, stick to) the very bottom of the wave potential 
troughs, so f i ^ := — ea/fc — e(p. Hence, 



The assumption that the number of trapped particles is 
fixed implies that (i) those are trapped in autoresonance 
deeply^ such that they do not become untrapped when the 
wave parameters evolve; (ii) also, it is implied that there 
are no passing particles close to the resonance, so that no 
additional trapping can result from the wave evolution. 
(As models, corresponding distributions already proved 
useful for understanding paradigmatic effects driven by 
trapped particles USHU]; yet, they can also form natu- 
rally as waves evolve [IS].) Then, assuming also that the 
wave envelope is smooth and evolves slowly enough, an 
adiabatic Lagrangian density £ in the GO approximation 
can be inferred from the general formalism reported in 
Paper I. For simplicity, we will focus on ID nondissipa- 
tive electrostatic waves in nonmagnetized plasma here. 
In this case, £ reads as (Paper I) 



o = 






a 
16^ 



-^ns(Hs)/,, 



(2) 



which describes both the quasistatic potential Cp (if any) 
and also the wave, to be characterized by the amplitude 
a, the frequency w, and the wave number k. The summa- 
tion in Ec^. ([2| is taken over different species; n^ are the 
corresponding densities averaged locally over the wave os- 
cillations; "Hs are the OC energies, and the angular brack- 
ets denote averaging over the particle distributions fs- 

We assume that oj is large enough, such that no ions 
are trapped and that the ion quiver motion is insignifi- 
cant. Then, 



ni{'Hi)f, = riiUi = S)i -I- ernip, 



(3) 



where S^i is the ion thermal energy density, e > is 
minus the electron charge, and n^ is the ion density, which 
is equal to the electron unperturbed density no. (For 
clarity, we take the ion charge to be e.) Suppose also that 
the number of trapped electrons is small enough so that 
the wave can be considered monochromatic; for specific 
conditions see Refs. 20, 26]|2J. Then, provided that all 
passing electrons are nonresonant, thus undergoing linear 
oscillations, one can write (Paper I) 



n'i'Hu^i'^)^ = ^e - (e - 1) ^ - en(^)^, (4) 



a 
16"7r 



where S^^ is the electron thermal energy density, and 
e(a;, fc) is the linear dielectric function. Finally, trapped 



(•Hi*))/, = -eni*)a/fc - en'^^Cp - m^n^^u^/2. (6) 



Combining the above equations yields 



ea^ iruTui'^ 

2. = h eaa + 

Idtt 



2fc 

Stt 



;[no 



rSP^]^ + eak^. (7) 



Here we omitted the insignificant thermal energies, 
dropped the index e, and introduced a = rv*' /k, which is 
proportional to the number of trapped electrons within 
the local wavelength. In this paper, we will assume 
a{x,t) = const for simpHcity; otherwise, see Paper I. 

In what follows, we will also assume that the wave is 
close to stationary at each x. Then, to prevent the qua- 
sistatic field build-up, the total electron current must re- 
main approximately zero (assuming that the ion current 
is zero). Hence, the fiow velocity of passing electrons, Vb; 
is estimated as [T8l [19] 



K) 



iW/n(P). 



(8) 



We will assume that such Vb does not affect the plasma 
dispersion, the condition being that the right-hand side 
of Eq. ((sl) remain small compared to other characteristic 
velocities in the system. In particular, we will require 
that it be small compared to the linear group velocity 
VgQ. (Remember that VgQ <^ vt ^ u, where vt is the 
passing-electron thermal velocity; see also Seclvl) Thus, 



N-: 



.(*) 



<1 



will be assumed, where 



K — kXz) <C 1, 



(9) 



(10) 



and Xd — vt/^p is the Dcbye length. Since the bulk 
plasma is supposed to be cold, one can then take [3S] 



e = 1 



w' 



a;2 — ik'^v^ ' 



(11) 



where Wp = [ATTn^^'e^ /mY''^ is the plasma frequency. 



III. WAVE EQUATIONS 

Varying £, with respect to the wave variables yields 
wave equations in the geometrical-optics approximation 
as discussed in Paper I (see also Appendix [A|. In partic- 
ular, SaS: = yields Eq. Q, or 



e{uj, k) + 87rcre/a = 0, 



(12) 



which defines the wave NDR, w = uj{k,a). Using 
Eq. (11), one obtains then (cf. Paper I) 



w^ = 



CO' 



1 + 2u;f/uj% 



+ Sk'^v^, 



(13) 



where uit — [iTTn'^^^e^/mY^^, and uje = {eak/mY^'^ 
is the characteristic frequency of the trapped-particle 
bounce oscillations in the wave troughs. Notice that the 
sinusoidal-wave approximation requires |20j 



■& = UJ^/UJ% < 1. 



(14) 



Thus, even with the nonlinear corrections, within our 
model one still has the approximate equality w^ « w^; 
i.e., studied here are weakly nonlinear Langmuir waves. 
Further, notice that k = dx£, and w — —dt£,, where ^ is 
the wave phase. Then 3^2 = gives the ACT, 



dtl + d,J = 0, 
where X = £,^ is the action density, 

9 



Wtt 



and J' = — £fc is the action flux density. 



J 



£ka 
' Wtt 



eif 



(15) 



(16) 



(17) 



As argued in Appendix [^ the contribution of eip in 
Eq. (171 can be neglected in the paradigmatic regimes 



that we consider here. (For effects of prescribed nonzero 
(f, which also flow from Eq. ( 15 1, see Ref. [29 .) This elim- 
inates (f from the wave equations altogether, so the wave 
Lagrangian can henceforth be used in a reduced form, 



£(a, oj, k) — e(w, k) — h ecra + 



16tt 



2k 



yielding, in particular, that 



a 

Siruj 



mau, J ; 



Skv^a'^ 

87rcj2 



mau 



(18) 



(19) 



Remarkably, the contribution of autoresonant trapped 
particles to £ is of lower order in a than the leading 
(first) term. Compared to common nonlinear waves, this 
is quite unusual. Moreover, the third term in £, which 
is due to the trapped-particle OC inertia, is completely 
independent of the field amplitude and yet cannot be 







(a) v= 1x10"'' 

(b) v=3xlO"^ 

(c) v/ = 5xlO"^ 



n/tic 



FIG. 1: ID stationary wave with trapped particles in inho- 
mogeneous plasma. Shown is the normalized amplitude of 
the wave potential, <j)/<j)c, vs. the plasma normalized density 
n/ric = Op for different u. Here (f>c = 20, with Q = 0.1 taken 
as an example; ric is the critical density. Dashed is the linear 
solution {u = 0) and the location of the linear cutoff {n — Uc). 



dropped as it still depends on w and k. Below, we show 
how certain paradigmatic effects result from this inertia 
term, specific to waves with trapped particles. Notice 
also that, at least qualitatively, those effects are not lim- 
ited to the model of deeply trapped particles that we 
adopt. This is because the a-independent term in £ orig- 
inates from the second term in Eq. ([5]) , which is indepen- 
dent of £e and thus yields a contribution invariant to 
the trapped-particle distribution. In particular, the ap- 
plication of the Whitham's approach in Ref. [T3] is hence 
invalidated I3QJ. and some other existing theories must 
be revised, as we will argue in Sec. |VI| 



IV. 



WAVE TRANSFORMATIONS IN PLASMA 
WITH VARYING PARAMETERS 



Now that £ is known explicitly, one can apply the 
ACT to infer the wave evolution in plasma with parame- 
ters varying in space and time. In particular, integrating 
Eq. ( 15 ) over the volume predicts the conservation of the 
wave total "number of quanta", JXdx = const. Com- 
bined with the NDR, for a homogeneous wave this yields, 
for instance, a = a{uj{t)) in the same manner as for linear 
waves discussed in Refs. [551 [3T1 [32] . Below, we consider 
another example, namely, stationary wave propagation 
in inhomogeneous plasma. Unlike in Refs. [TTHT^ , where 
an ad hoc solution of the Vlasov equation was employed 
to address a similar problem, we will show that our ACT 
yields the same results straightforwardly. 

For a stationary wave Eq. (15) gives J^ — const. Since 
one also has dxUJ = in that case [Eq. (B2)], this is un- 
derstood as conservation of the wave energy flux 11 — loJ 
(cf. Appendix [B]), or 



n(p) -I- n(*) = const, 



(20) 



where we introduced the densities of energy fluxes carried 



by passing particles and trapped particles, respectively: 



(cf. also Refs. [H] 



tl\, UrpU 

Sttlo 



n(*) = - n(*)mu3 



5]). If the ratio 



(21) 



(22) 



is small, one can anticipate that the wave is in the linear 
regime and is not affected significantly by trapped parti- 
cles. However, S can become large as the plasma density 
varies along the ray trajectory. Then the wave dynamics 
becomes essentially nonlinear, which is seen as follows. 

Let us introduce Qp{x) — ujp{x)/uj and x = kXx ~ k, 
where Xt = vx/uj ~ Xjj; also, v = a /{iiIcXt) ~ A/", with 
Tic being the critical density, and (f) — ea/{mjjOVT)- Then 



Eq. (20) rewrites as 



3>f(/)2 + v/>? = 39, 



(23) 



where O is some constant determined by boundary con- 



ditions. Also, we can substitute x from Eq. ( 13 1, namely, 
3^2 

Then the equation for (j) reads as 



{l-nl) + 2vl<t>. 



(24) 



V3 



1 



ni 



2vl(j)- 



1 



^l 



2v/(t) 



= e. 



(25) 



Hence, the effect of trapped particles is determined by 
two parameters: S", defined [in agreement with Eq. (22|] 
as the ratio of the terms on the left-hand side here, and 
also cr, which is a measure of how x is affected by the 
nonlinearity [cf. Eq. (24|]; namely. 



v/x^ 






(26) 



Suppose that a wave is launched from a subcritical 
region {Vtp < 1) and initially is close to linear (5* ^ 1, 
^ <C 1), the boundary condition thus being w xq(J)q. 
Then, at first, (/) ~ 4'oi^o/ ^)^ , so i;^ remains constant. 
Yet S grows as >c~'^, and eventually one gets S ^ 1, at 
som e X ~ Xif !34| . Beyond that point, the second term in 
Eq. (23) dominates, resulting in constant x^ « zy/(30). 
Hence Eq. E^ yields ^1^1 + v/Q = 2v/(j), or 



2v 



9? 

p 



1 + i^/e' 



(27) 



As flp grows further, (j) starts to decrease then (Fig.nl). 
Finally, the wave reaches x = Xc, where 51^ = 1 [and thus 
(;{xc) ~ 1]. At that location, (p equals about 26 and then 
continues analytically into the overcritical region. 

Hence, in contrast with a linear wave (for which a GO 
solution with nonzero H would be impossible there), a 
nonlinear wave loaded with trapped particles can, in prin- 
ciple, penetrate overcritical plasma, in agreement with 





FIG. 2: (Color online) Schematic of the parameter domain as- 
sumed for Sec.^ (a) in space («:, ■&, Qe), (b) in space (k, i9, S). 
Combined here are the following assumptions: the plasma is 
cold [Eq. ( 10 1] , the wave is sinusoidal [Eq. ( 14 \] and weak 



enough [Eq. (A7l], the quasistatic field due to trapped par- 
ticles is negligible [Eq. (Alll]; also, the bulk motion and the 



nonlinear effects are w eak , i.e., Vb ^ Awj, <^ Vgo [Eq. (30 1; 
see also Eqs. ^ and ([28}]. The inequalities ^ and ( |Al2 | 
reduce to Vo ^ Vgo and thus are satisfied automatically. 



Refs. p7H19j . In particular, since the inhomogeneity 
scale does not enter the above equations explicitly, it can 
always be chosen large enough, so as to ensure the valid- 
ity of the GO profile we derived. However, notice that 
such a wave will propagate adiabatically for limited time 
only. This is because S{x > Xc) > (i^9)^^/^ >• 1, in 
which case the wave is unstable, as we will now discuss. 



V. PULSE PROPAGATION 

Let us consider the envelope dynamics, assuming, for 
simplicity, that the plasma is homogeneous and station- 
ary. Since the ACT [Eq. (15)], serving as the envelope 
equation, contains only first-order derivatives of a, the 
diffraction is hereby neglected. Within this approxima- 
tion, a linear pulse would conserve its shape, traveling 
at the linear group velocity VgQ. A nonlinear wave, in 
contrast, will undergo distortion, particularly because it 
can have two group velocities Vg. 



A. Group velocity splitting 

Specifically, following Refs. [Tl[2], we define Vg here as 
the velocities of information (rather than that of the en- 
ergy, like in Refs. [55H55] ). which propagates along char- 
acteristics of the wave equation (Appe ndix [B |). Hence, Vg 
are found explicitly from Eqs. ( B21|)- ( [B24[ ), via substi- 
tuting Eq. pSl) for £. Using Eqs. (Iiop and pll, after a 



tedious yet straightforward calculation, one gets then 



VgO 



1 + S^±g 
1 + 3S-dK^ ' 



(28) 



where Vgo w 3kvt, and 

g = ^eVS {1/2 ~S), 



(29) 



with riE = loe/'^p', in particular, S — ^/(3ri|;K^). We as- 
sume that the nonhnear corrections to VgQ are smaU, i.e., 

Sd-^1, <7 < 1, (30) 

where the former inequahty ensures that the denominator 



in Eq. ( 28 1 is also close to unity. (For a summary of our 
assumptions, see Fig. [2J) Hence, 



VgO 



i±nEVsii/2-s) 



(31) 



with a characteristic shape of Vg{K) shown in Fig. ^ 

Although small, the second term in Eq. ( [3l| ) is re- 
tained because it is responsible for essentially nonlinear 
effects. In particular, consider the case when S < 1/2. 
Then g is real, causing the group velocity splitting by 
Avg — 2vgQg. This means, for example, that a gen- 
eral modulation imposed on the wave profile eventually 
splits into two signals propagating with different veloci- 
ties [Fig. Wa)]. Each signal may then evolve further, if 
having a finite spread of a and thus of Vg too; however, 
such a signal will be comprised of characteristics that all 
correspond to the same sign in Eq. (JSTJ) , so further split- 



ting per se will not occur. As we remind in Appendix [Bj 
the pulse splitting is an inherent feature of all nonlinear 
waves, as well known in classical hydrodynamics [21 [39] 
and also observed in plasma physics experiments [5J I49| . 



B. TPMI 



In contrast, at S > 1/2, Eq. ([31]) yields no real solu- 
tions for Vg. Then, the wave is TPMI- unstable, as illus- 
trated in Fig.llFb). One can also assess the wave stability 
using another approach, which yields the TPMI rate ex- 
plicitly. Instead of searching for characteristics, take 

a = flo + Aa, Lu = ujQ + Acj, k = ko + Afc (32) 





FIG. 3; (a) Nonlinear group velocities, Vg and Vg [Eq. (31l] 



vs. K, for sample Q,e and -d; the dashed line shows the linear 
group velocity Vgo- At S > 1/2, corresponding to «: < ks = 
r2^^\/2i9/3, no real solutions exist for Vg, rendering the wave 
unstable, (b) Close-up at k ~ us- 
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FIG. 4: (Color online) Evolution of the perturbation A(ao) = 
— O.OSflo exp(— a;'*/^'*) to a homogeneous wave with the initial 
amplitude ao. The solution is obtained by numerical integra- 
tion of Eqs. ([B1])-(1B3|, with £ taken from Eq. ([l8|, for the 
same parameters as in Fig. ^and i = 20Ao. Shown is A(a'^) 
(arbitrary color scaling), vs. t and x in the frame moving 
with the linear group velocity Vgo; the units are a;j7^ and Xr>, 
correspondingly, (a) n — 0.2, so 5* = 1/4; the wave is stable, 
resulting in signal splitting, (b) k = 0.1, so 5* = 1; the wave 
is TPMI-unstable. 



and assume that the perturbations, denoted by A, are 
small compared to the corresponding homogeneous pa- 
rameters, denoted by 0. Then, one can linearize the GO 



equations [Eqs. (B1)-(B3|] and find the dispersion rela- 
tion Aa;(AA;), which is linear in the sense that it is inde- 
pendent of Aa (albeit not of oq). Specifically, under the 
same assumptions as for Eq. (31 1, one can show that 



ALu^Akvga{l±g) 



(33) 



At S' < 1/2, one has real Aw, and the signal velocity 
J that we introduced 
linear group velocity 



Vg that we introduced earlier is recovered as the effective 



djAuj) 
d{Ak) 



VgO 



:i±5)- 



(34) 



At S > 1/2, one has two w complex conjugate to each 
other. Then the wave is unstable, as predicted earlier; 
yet now we get a formula for the TPMI rate too, 

7 « Afc Vgo^E ^SiS-1/2). (35) 



VI. DISCUSSION 

A. Stability criterion 

The TPMI threshold also can be inferred from the gen- 
eralized Lighthill's criterion (GLC), which states that a 
wave is stable when ujx{k,I)J'^.{k,I) > (Appendix [B| . 
To see this, let us first rewrite X as 



I: 



maujr, 



SiTiOr, 



(36) 



from where one gets a = a{k,I), with ax > 0. Since 
iOa{k,a) > 0, the latter yields uJx{k,X) > 0, so the GLC 
takes the form Ji^ik^I) > 0. Further, using that 



J 



Sttlo^ 



2fc2 



we obtain 



Jk 



34a^ 

47rw2 



kak \ 
a I 



(37) 



(38) 



Since kuk/a « 3Sk^ < S [from Eq. (36l], the GLC 



hereby rewrites as 5 < 1/2, in agreement with the crite- 
rion that we found earlier from Eq. ( |31[ ). 

Notice that waves with other distributions of trapped 
electrons also can be studied similarly. (Of course, a still 
must be conserved, which is not the case, e.g., for dis- 
tributions nonzero near the trapping boundary; cf. Pa- 
per II.) For example, if Wa(fc, a) < 0, the GLC would read 
a.s S > 1/2, assuming that Eqs. (36) and (37 1 still hold 
approximately. However, remember that the GLC, and 
the underlying Whitham's formalism overall, neglects 
diffraction and assesses the stability of only adiabatic per- 
turbations with small enough Afc. Even when those are 
stable, at larger Ak a more general SI [2TI - f23l |40] may 
develop [50], for which the TPMI, when present, can be 
considered as the adiabatic limit. In particular, notice 
that it is the general SI rather than the TPMI that is 
responsible for the effects reported in Ref. [40] . 



VII. CONCLUSIONS 

In this paper, the evolution of adiabatic waves with 
autoresonant trapped particles is analyzed under the as- 
sumption that the action distribution of these particles is 
conserved, and, in particular, that their number within 
each wavelength is a fixed independent parameter of the 
problem. Particularly, ID nonlinear Langmuir waves 
with deeply trapped autoresonant electrons are addressed 
as a paradigmatic example, illustrating the qualitative 
physics that is also expected for other distributions (ex- 
cept when dynamic trapping and detrapping become es- 
sential). For a stationary wave, tunneling into overcriti- 
cal plasma is explained from the standpoint of the action 
conservation theorem. For a nonstationary wave, quali- 
tatively different regimes are realized depending on the 
initial parameter S, which is the ratio of the energy flux 
carried by trapped particles to that carried by passing 
particles. At 5* < 1/2, a wave is stable and exhibits group 
velocity splitting. At S > 1/2, the trapped-particle mod- 
ulational instability develops, in contrast with the ex- 
isting theories of the TPMI yet in agreement with the 
general SI theory. Remarkably, these effects are not cap- 
tured by the nonlinear Schrodinger equation, which is 
traditionally considered as a universal model of wave self- 
action but misses the trapped-particle OC inertia. 
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B. Comparison with the existing theories 

The SI is treated systematically in Ref. [20] from where 
a rate can be inferred that matches our Eq. ( 35 ) under 
the appropriate conditions (AppendixO. Yet, our results 
are in drastic variation with the traditional |41j models of 
the TPMI [9HTT] , for those predict that a wave is always 
stable when resonant electrons remain deeply trapped [S] ■ 
The discrepancy is due to the fact that Refs. [SHU] rely 
on the NLSE, which neglects the contribution of the a- 
independent term in £ and thus is generally inapplica- 
ble to waves with trapped particles [32]. A more de- 
tailed comparison would require accounting for diffrac- 
tion (which is out of the scope of our model), but it is 
already seen that the NLSE-based models do not apply 
at large S. (Interestingly, disagreement between those 
models and numerical simulations has already been re- 
ported before; see, e.g., Refs. [TT1I43].) 

In connection to this, notice that at k > 0.2, which is 
usually associated with large enough n^*) [5S], S > 1/2 is 
achieved easily, not to mention that SI may develop too. 
Thus, the effect of the trapped-particle OC inertia on 7 
that we discuss here is important for assessing the wave 
stability under practical conditions. Also notice that the 
actual NDR may be nonlocal (Paper II) , contrary to the 
tacit assumption used commonly. Hence, further inves- 
tigation may be needed to adequately model the TPMI 
and related effects. 
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Appendix A: Electrostatic potential 

Let us discuss the specific conditions under which the 
electrostatic potential <f can be neglected in Eq. (17 1. 
First of all, the equation for (p is obtained from (5<^£ = 
[Eq. ([7])] and represents the usual Poisson's equation. 



dl 



,f 



-Aire [riQ 



,(p) 



,(*) 



(Al) 



where uq is determined by initial conditions, and n^*' is 
"frozen" into the wave, due to tr = const. To find n^P\ 
consider the momentum equation for passing electrons, 

m (dtVo + Vo d,Vo) = d^{e^ - $) - T9, lnn^P\ (A2) 

where we assumed, for simplicity, that the plasma is 
isothermal and included the ponderomotive potential, 
$ « e^a-^/{4mujp) (in contrast with Ref. [TS], where <& 
was left out). 

Since ^/uiVq ^ i?^^ ^ 1, we can neglect the convec- 
tive term on the left-hand side immediately. The other 
term can be estimated from 



mdtVo ^ mVgodxVo, 



(A3) 



with tilde henceforth denoting perturbations to station- 
ary values. From the continuity equation, 

a^n^P) + a^n'^Vo] = 0, (A4) 

one gets that AVg ^ Vgon^P^/n'-P^; hence, 

mdtVo - mv^go d^fi^P'^/n'^P^ <^T d.^fi^P^ /n^P\ (A5) 

This allows one to neglect the inertia term compared to 
the pressure term in Eq. ( A2 ) , thus yielding 



d^[-e(p + ^ + T\nn'^P^] =0. 
For weak enough waves, with $ ^ T, or 

one can thereby estimate the maximum (p as 
e(p "^T ^ mu . 



(A6) 



(A7) 



(A8) 



For the purpose of Sec. IV Eq. ( |A8[ ) is already enough 
to neglect e(p in Eq. (17). For the other regimes that we 



disc uss, one may assume e(p <^ T, yet what matters in 
Eq. ( |17| now are only variations of mu^ , which may not 
be smaller than e(p; thus, Eq. (A8| becomes insufficient. 
In this case, let us revert to Eq. ( Al ) and substitute 



,(p) 



^(rtg(e^-*)/T ^ ^(P) ^^ ^ (g - _ ^yj.^ (Ag) 



from Eq. ( A6 1 ; hence 



Xldl,{e^)^e^~'P + Tn^''>/nlJ 



(p) 



(AlO) 



Ap) 



where tiq is the unperturbed density of passing elec- 



trons. The latter term in Eq. (AlO) is negligible, due to 



$ „(p) 



< Sk" < 1, 



(All) 



provided that S is moderate. (Remember that we are 



mostly interested in regimes with S ^ 1.) At A^, d^^ <C 1, 



one gets then eip ~ <&. On the other hand, 

'AA< 1 



Cfca^/IGTT 



(A12) 



[Eq. i9l]. Thus, again, e(p is negligible in Eq. (17) 



Appendix B: Nonlinear group velocities and wave 
stability 

In this appendix, we concisely restate, to avoid am- 
biguity, the concept of the group velocity for linear and 
nonlinear waves, as introduced originally by Whitham 
[B H] (see also Refs. [351 IUHlZj), and discuss how it is 
connected with the wave stability. 



1. Conservation laws and flow velocities 

For simplicity, we will assume a ID system, the gen- 
eralization to the case of multiple dimensions being 
straightforward. In the GO approximation, the wave is 
completely described by its amplitude a, frequency w, 
and the wave number fc, so the Lagrangian density can 



be taken in the form £(a, uj, k; t, x). Then, like in Sec. Ill 
the wave equations read as 



£a=0, 


(Bl) 


dtk + dxUj = 0, 


(B2) 


dt£^ - d^£k = 0. 


(B3) 



Equation (Bl) defines the NDR, u — Lu{k,a). Equation 



(B2) is the consistency condition, satisfied due to 
w = -dt£., k = (9a;^, 



(B4) 



whe re ^ is the field phase. [Notice, in particular, that 
Eq. ( B2 1 can be understood as the continuity equation for 
wave crests, with k being the crest density, and u — ku 



being the crest flux density.] Finally, Eq. (B3) is the 
ACT, also known in the form ( 15 ), or 



dtT + d^{viI)^Q, 



(B5) 



with vx — —Zkl^uj serving as the action flow velocity. 

Suppose now that the medium is stationary, i.e., £ does 
not depend on t exp licit ly. Then, 9t£ = £^ dtto -f £fe dtk, 
where we used Eq. (Bl). Thus, 



5t(w£^-£)=cjat£„- 

= LO dt&to - 
= w9a;£fe - 

= ijjdx^k - 



£„ dtuj - 
2k dtk 
£fc dtk 

■ £fe dxU! 



&UJ dtio - Zkdtk 



(B6) 



where Eqs. (B2) and (B3| were employed. Similarly, if £ 



does not depend on x explicitly, one gets 
dx{kU-2)=dt{kZ^). 



(B7) 



Equations (B6) and (B7| represent the conservation laws 



for the wave energy and momentum and can be written as 
dts + dx{vee) = 0, dtp + dx{vpp) = 0, (B8) 

where we introduced 

£ = uj£,^ - £, p = fc£^ (B9) 

for the energy density and the momentum density and 
w£fc A:£fe - £ 



c^£^-£' 



fc£. 



(BIO) 



for the corresponding flow velocities. In particular, notice 
that the energy flux density II = v^e and the momentum 
flux density V = Vpp are then given by 



n 



-LuHk, V^Z-kZk 



(Bll) 



In general, vx, v^, and Vp are all different from each 
other. The exception is the linear regime, which is de- 
fined as the regime when uj{k), inferred from Eq. (Bl|, 
is independent of a. The latter is possible if £,a has the 
form £,a = 23(w,A:)Ao, where A is some function such 



one can rewrite the ACT as 



that Aa is nonzero; hence £ — S)(w, k)A. Since Eq. (Bl) 
thereby reads as S(i^, k) = 0, one has £ = 0, so 



vi 



-£fe/£. 



(B12) 



From differentiating £(a, w(A:), fc) = with respect to fc, 
one gets —£,k/S.ui = ^k = Vgo, with the latter known as 
the linear group velocity. Therefore, in the linear regime. 



VI 



"so- 



(B13) 



Notice also that a pulse is usually linear at its front and 
tail (except in the presence of trapped particles), since 
the field is weak there. Hence, it is only within the pulse 
that Eq. (B13) can be violated, due to nonhnear effects. 



2. Nonlinear group velocity 



By analogy with the linear case [Eq. (B13)], the non- 
linear group velocity Vg is often defined as «£ too [55H55] , 

or as Wfe with the derivative taken at fixed £/w, since 



{'^k)s./u 



dki^H_ 



£,u - S./U 



(B14) 



However, this generalization is arbitrary, and other def- 
initions, such as Vg = Vp or Vg = vi, would be equally 
justified. [In fact, the latter would be more fundamental, 
because the ACT holds also in nonstationary medium, 
unlike the energy conservation law.] More consistently, 
Vg is defined as the velocity of information, i.e., the ve- 
locity on characteristics [3 [5] . Below, we restate how it 
is calculated in the general nonlinear problem. 



To find characteristics of Eqs. (Blj)-(B3), let us con- 
sider traveling-wave solutions, witi 
locity Vg yet to be found. 



;h the propagation ve- 
in other words, let us search 
for solutions in the form where all the wave variables are 
expressed through a single variable CC^;,^ ^ x — X{t), 
such that dtX — Vg] then, d^ — d,^ an d dt — —Vgd^. In 
particular, one thereby gets from Eq. (|B2|) that 



— uj' /k' = (ifcOj, 



(B15) 



where primes denote d^, and dk is taken in the sense 
that u! — uj{k, a{k)). [Notice that the latter local relation 
holds on a characteristic only, and generally there are two 
branches of a;(fc, a{k)) corresponding to the two different 
types of characteristics that we will find.] Hence, under- 
standing the nonlinear group velocity as the characteris- 
tic velocity represents a natural generalization of VgQ. 
To actually find Vg, we proceed as follows. Using that 

dt£,^ = £u;a dta + £^^ dtuj + 2^k dtk, (B16) 

dx£-k = £fea d^a + £kk dxk + £fccj d^uj, (B17) 



Wgl-Cojaa' 



£fcaa' + £fefefc' + £fc^w' = 0. (B18) 



Here a' can be derived from Eq. (Bl ), after differentiating 
the latter with respect to C,: 



= d^Za = 2,aa a + S^aui^' + &akk' ■ 

Specifically, one gets 



(B19) 



da Zau^Vg + £afc 


(B20) 


so Eq. (|B18| rewrites as follows 




pv] + 2rvg + q = 0, 


(B21) 


where we introduced 




P ^ ZaaZu:uj ~ ^oja; 


(B22) 


r =^ ^aa^uk '^Lja^ka^ 


(B23) 


q = Zkk^aa - £fca- 


(B24) 



Since Eq. ( B21 1 is a quadratic equation for Vg, there are 



generally two group velocities different from each other 
(regardless of the type of nonlinearity) , causing signal 
splitting. The exception is the linear regime. In that 
case, it is convenient to use A instead of a (Sec. Bl|; 



then the same equations hold, if da is replaced with Oa- 
On the other hand, £,aa — 0, so one obtains 

v^ZIa + 2vg£^AQkA + -CL - (B25) 

(cf. also Ref. [35] )i yielding that the two roots coincide: 

Vg = -£fcA/£^A - -£fc/£^ = vgo. (B26) 

Thus, we again see that the nonlinear Vg, defined as the 
characteristic velocity, equals VgQ in the linear limit. 

3. Stability criterion 

If r^ < pq, there are no real solutions for Vg, so no 
stable envelope is possible in this regime. This means 
that amplitude modulations will grow with time, i.e., the 
wave is modulationally unstable, unless r^ > pq. Below, 
we will put this criterion in a yet different form 2 , which 
we will need in the main text. 

Let us choose the action density Z = £1^ to serve as an 
independent variable instead of a. Also, using the NDR, 
exclude w from the list of independent variables; hence, 

ij^uj{k,I), a^a{k,I), J = J{k,I), (B27) 

where J^ = — £fc is the action flux density. In particular, 
notice that £ takes the following form: 



il{a{k,I),uj{k,I),k) = A{k,I), 



(B28) 



so Afc = lijjk — J and Ax = Xwx, due to Eq. (Bl ). Now 



Eqs. (B2) and (B3) on characteristics read as 



-Vgk' + Wfcfc' + cjil' = 0, -Vgl' + Jkk' + Jil' == 0. 
From Afei = Axfc, it follows that w/j — Jx-, so one gets 



Vg=UJ],± \/uJxJk- 



(B29) 



(In a linear wave, ix>x — 0, so there is only one group 
velocity, Vg — VgQ.) Hence, the wave is stable if 



uJxJk > 0. 



(B30) 



In particular, when wj is small, one can substitute in 
Eq. (B30I the lowest-order approximation for J', yield- 
ing LOxJ^k = '^x^dkVgQ. Suppose that, approximately, 
I ex a^ (corresponding to the most common choice of 
a), so ujxl = ujaa/2. Then, required for stability is the 
condition coa dkVgo > 0. The latter is equivalent to that 
flowing from, e.g., Eq. (4) in Ref. [5] or Eq. (7) in Ref. [15] 
and also agrees with the Lighthill's original criterion for 
weakly nonlinear waves |46| . Hence, in the main text, 
Eq. (B30) is called the generalized Lighthill's criterion. 



Appendix C: TPMI as the adiabatic limit of SI 

Let us compare the TPMI rate that we found in Sec.[V| 
with that flowing from Ref. [5D] for the more general SI 
[2Tti23l l40]. Specifically, we will consider the hmit 5* > 1, 



where Eq. (23) of Ref. [50] applies, reading as 
A (SAT)^ 



1- 



Alj^ [ACu + Ak{l-TW 



= 0. 



(CI) 



Here A =_51|, T = 3^2, Alu = {Acu - Akuo)/ujo, uq = 
LUo/ko, Afc = Ak/ko, and S is defined as in our paper. 

Equation ( CI ) describes four eigenmodes, two of which 
correspond to oscillations at Aw sa ±A^'^ and two others 
have yet lower frequencies. We will assume that the for- 
mer are of zero amplitudes (remember that our theory 
applies only at time scales large compared to ui]^ ; see 
Paper I) , so those are the lower- frequency modes that we 
will consider. Hence we take A<D^ <C A, in which case 



the first term in Eq. (CI) can be neglected, yielding 



[Aw + Afc(i-r)]' = -- 



A 



[SATy 



(C2) 



Treating the right-hand side as a perturbation, one gets 
Au K Ak{T - 1) ± iAkSrVA (C3) 

(here we also used that T ^ 1), which is equivalent to 

Aw w Ak vgo ± i{Ak/k) ujpSTVA. (C4) 

Since TA^^^ — SflEU^, the instability rate 7 then equals 
j^AkvgoQES, (C5) 

which precisely matches our Eq. ([35]) taken at 5 » 1. 
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